Formation of large-scale semi-organized structures in turbulent convection 
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A new mean-field theory of turbulent convection is developed by considering only the small-scale 
part of spectra as "turbulence" and the large-scale part, as a "mean flow", which includes both 
regular and semi-organized motions. The developed theory predicts the convective wind instability 
in a shear-free turbulent convection. This instability causes formation of large-scale semi-organized 
fluid motions in the form of cells or rolls. Spatial characteristics of these motions, such as the 
minimum size of the growing perturbations and the size of perturbations with the maximum growth 
rate, are determined. This study predicts also the existence of the convective shear instability in 
a sheared turbulent convection. This instability causes generation of convective shear waves which 
have a nonzero hydrodynamic helicity. Increase of shear promotes excitation of the convective shear 
instability. Applications of the obtained results to the atmospheric turbulent convection and the 
laboratory experiments on turbulent convection are discussed. This theory can be applied also for 
the describing a mesogranular turbulent convection in astrophysics. 

PACS numbers: 47.65.+a; 47.27.-i 



I. INTRODUCTION 



In the last decades it has been recognized that the 
very high Rayleigh number convective boundary layer 
(CBL) has more complex nature than might be reck- 
oned. Besides the fully organized component naturally 
considered as the mean flow and the chaotic small- 
scale turbulent fluctuations, one more type of motion 
has been discovered, namely, long-lived large-scale struc- 
tures, which are neither turbulent nor deterministic (see, 
e.g., 0, | | I 1,0,0,00,0). These 

semi-organized structures considerably enhance the ver- 
tical transports and render them essentially non-local in 
nature. In the atmospheric shear-free convection, the 
structures represent three-dimensional Benard-type cells 
composed of narrow uprising plumes and wide down- 
draughts. They embrace the entire convective bound- 
ary layer (~ 2 km in height) and include pronounced 
large-scale (~ 5 km in diameter) convergence flow pat- 
terns close to the surface (see, e.g., |^], and references 
therein). In sheared convection, the structures represent 
CBL-scale rolls stretched along the mean wind. Life- 
times of the semi-organized structures are much larger 
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than the turbulent time scales. Thus, these structures 
can be treated as comparatively stable, quasi-stationary 
motions, playing the same role with respect to small-scale 
turbulence as the mean flow. 

In a laboratory turbulent convection several organized 
features of motion, such as plumes, jets, and the large- 
scale circulation, are known to exist. The experimen- 
tally observed large-scale circulation in the closed box 
with a heated bottom wall (the Rayleigh-Benard ap- 
paratus) is often called the "mean wind" (see, e.g., 
||, [U| H H fH HI H, and references therein). 
There are several unsolved theoretical questions concern- 
ing these flows, e.g., how do they arise, and what are their 
characteristics and dynamics. 

In spite of a number of studies, the nature of large-scale 
semi-organized structures is poorly understood. The 
Rayleigh numbers, Ra, based on the molecular transport 
coefficients are very large (of the order of 10 11 — 10 13 ). 
This corresponds to fully developed turbulent convec- 
tion in atmospheric and laboratory flows. At the same 
time the effective Rayleigh numbers, KsS cS \ based on 
the turbulent transport coefficients (the turbulent viscos- 
ity and turbulent diffusivity) are not high, e.g., Ra*- 05 -* ~ 
Ra/(RePe), where Re and Pe are the Reynolds and Peclet 
numbers, respectively. They are less than the critical 
Rayleigh numbers required for the excitation of large- 
scale convection. Hence the emergence of large-scale con- 
vective flows (which are observed in the atmospheric and 
laboratory flows) seems puzzling. 

The main goal of this study is to suggest a mechanism 
for excitation of large-scale circulations (large-scale con- 
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vection). In particular, in the present paper we develop 
a new mean-field theory of turbulent convection by con- 
sidering only the small-scale part of spectra as "turbu- 
lence" and the large-scale part, as a "mean flow", which 
includes both, regular and semi-organized motions. We 
found a convective wind instability in a shear-free turbu- 
lent convection which results in formation of large-scale 
semi-organized fluid motions in the form of cells or rolls 
(convective wind). We determined the spatial charac- 
teristics of these motions, such as the minimum size of 
the growing perturbations and the size of perturbations 
with the maximum growth rate. In addition, we stud- 
ied a convective shear instability in a sheared turbulent 
convection which causes a generation of convective shear 
waves. We analyzed the relevance of the obtained results 
to the turbulent convection in the atmosphere and the 
laboratory experiments. 

Traditional theoretical models of the boundary-layer 
turbulence, such as the Kolmogorov-type closures and 
similarity theories (e.g., the Monin-Obukhov surface- 
layer similarity theory) imply two assumptions: (i) Tur- 
bulent flows can be decomposed into two components 
of principally different nature: fully organized (mean- 
flow) and fully turbulent flows, (ii) Turbulent fluxes are 
uniquely determined by the local mean gradients. For 
example, the turbulent flux of entropy is given by 



( SU ) 



-kt'VS 



(1) 



(see, e.g., |23|), where kt is the turbulent thermal con- 
ductivity, S is the mean entropy, u and s are fluctuations 
of the velocity and entropy. 

However, the mean velocity gradients can affect the 
turbulent flux of entropy. The reason is that additional 
essentially non-isotropic velocity fluctuations can be gen- 
erated by tangling of the mean-velocity gradients with 
the Kolmogorov-type turbulence. The source of energy 
of this "tangling turbulence" is the energy of the Kol- 
mogorov turbulence. 

In the present paper we showed that the tangling tur- 
bulence can cause formation of semi-organized structures 
due to excitation of large-scale instability. The tan- 
gling turbulence was introduced by Wheelon |24| and 
Batchelor et al. [^5| for a passive scalar and by Golit- 
syn |2(| and Moffatt (27]] for a passive vector (mag- 
netic field). Anisotropic fluctuations of a passive scalar 
(e.g., the number density of particles or temperature) 
are generated by tangling of gradients of the mean pas- 
sive scalar field with random velocity field. Similarly, 
anisotropic magnetic fluctuations are excited by tangling 
of the mean magnetic field with the velocity fluctuations. 
The Reynolds stresses in a turbulent flow with a mean 
velocity shear is another example of a tangling turbu- 
lence. Indeed, they are strongly anisotropic in the pres- 
ence of shear and have a steeper spectrum (oc fc -7 / 3 ) than 
a Kolmogorov turbulence (see, e.g., §8[ §|, © !!)• Thc 
anisotropic velocity fluctuations of tangling turbulence 
were studied first by Lumley p8| . 

This paper is organized as follows. In Section II we 



described the governing equations and the method of the 
derivations of the turbulent flux of entropy and Reynolds 
stresses. In Section III using the derived mean field equa- 
tions we studied the large-scale instability in a shear- 
free turbulent convection which causes formation of semi- 
organized fluid motions in the form of cells. In Section 
IV the instability in a sheared turbulent convection is 
investigated and formation of large-scale semi-organized 
rolls is described. Application of the obtained results for 
the analysis of observed semi-organized structures in the 
atmospheric turbulent convection is discussed in Section 
V. 



II. THE GOVERNING EQUATIONS AND THE 
METHOD OF THE DERIVATIONS 



Our goal is to study the tangling turbulence, in par- 
ticular, an effect of sheared large-scale motions on a de- 
veloped turbulent stratified convection. To this end we 
consider a fully developed turbulent convection in a strat- 
ified non-rotating fluid with large Rayleigh and Reynolds 
numbers. The governing equations read: 



0_ 

dt 



v • V v = -V 



-g5 + k(v), (2) 



i+vvW-vNi-iv.P^S), (3) 



where v is the fluid velocity with V ■ v = A ■ v, g is 
the acceleration of gravity, pQf v (y) is the viscous force, 
F K (S) is the heat flux that is associated with the molec- 
ular heat conductivity k, A = — p^ 1 V/?o is the density 
stratification scale, and N& = (^P )~ 1 VP — Pq 1 ^Po- 
The variables with the subscript "0" correspond to the 
hydrostatic equilibrium VPo = PoS, and To is the equi- 
librium fluid temperature, S — P/jPo — p/po are the de- 
viations of the entropy from the hydrostatic equilibrium 
value, P and p are the deviations of the fluid pressure 
and density from the hydrostatic equilibrium. Note that 
the variable S = 0/Oo, where O is the potential temper- 
ature which is used in atmospheric physics. The Brunt - 
Vaisala frequency, Oj, , is determined by the equation 
fll = — g ■ IS! i,. In order to derive Eq. (^) we used an iden- 
tity: -VP + gp = -po[V(P/po) + gS-PN b /po], where 
we assumed that \PN b /p \ < |g5| and |PN fc /p | < 
|V(P/po)|- This assumption corresponds to a nearly isen- 
tropic basic reference state when Nj, is very small. For 
the derivation of this identity we also used the equation 
for the hydrostatic equilibrium. Equations (|^) and (^|) are 
written in the Boussinesq approximation for V • v =/= 0. 



A. Mean field approach 

We use a mean field approach whereby the velocity, 
pressure and entropy are separated into the mean and 
fluctuating parts: v = U + u, P = P + p, and S = 
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S + s, the fluctuating parts have zero mean values, U = 
(v), P= (P) and S = (S). Averaging Eqs. @ and §) 
over an ensemble of fluctuations we obtain the mean-field 
equations: 



d_ 

dt ' 

d_ 

at 



U • V K 



U- V )S 



" V '(£) + ( A J - v j)<«<«3> 
-gS + f„(U), (4) 
-U-N 6 + (Ai- V0(*«i) 
f 

~7} 



V-F K (U,S) , 



(5) 



where pof^U) is the mean molecular viscous force, 
P«(U, S) is the mean heat flux that is associated with 
the molecular thermal conductivity. In order to derive 
a closed system of the mean-field equations we have to 
determine the mean-field dependencies of the Reynolds 
stresses /^(U, S) — (iij(t, x)itj(t, x)) and the flux of en- 
tropy $i(U, S) — (s(t,x)ui(t,x)). To this end we used 
equations for the fluctuations u(t, r) and s(t, r) which 
are obtained by subtracting equations (^) and (0) for the 
mean fields from the corresponding equations (Q) and (|^) 
for the total fields: 



du 

~dt 

ds 

dt 



-(U- V)u- (u- V)U- V( — 

-gs + U N , (6) 
-u-(N b + V5)-(U-V)s + S w , (7) 



where Un = ((u ■ V)u) — (u • V)u + f w (u) and Sn = 
((u-V)s) — (u- V)s— (1/To) V-F K (u, s) are the nonlinear 
terms which include the molecular dissipative terms. 



B. Method of derivations 

By means of Eqs. (^J) and (0) we determined the depen- 
dencies of the second moments /y (U, S) and 3>j(U, 5) on 
the mean-fields U and S. The procedure of the derivation 
is outlined in the following (for details see Appendix A). 

(a). Using Eqs. (0) and (0) we derived equations for 
the following second moments: 



/y(k) = L(ui,Uj) 
F(k) = L(s,w) , 
H(k) = L(s,s) , 



$i(k) = L(s,0 , (8) 
G(k)=L(w,w), (9) 
(10) 



where L(a,b) — (a(k)6(— k)), cl> = (Vxu) z , the accel- 
eration of gravity g is directed opposite to the z axis. 
Here we used a two-scale approach. This implies that we 
assumed that there exists a separation of scales, i.e., the 
maximum scale of turbulent motions Iq is much smaller 
then the characteristic scale of inhomogeneities of the 
mean fields. Our final results showed that this assump- 
tion is indeed valid. The equations for the second mo- 



ments (|)-@ are given by Eqs. ( |A15| ), JA16| ) and flA2l| )- 
( |A23| ) in Appendix A. In the derivation we assumed that 
the inverse density stratification scale A 2 <C k 2 . 

(b). The derived equations for the second moments 
contain the third moments, and a problem of closing the 
equations for the higher moments arises. Various ap- 
proximate methods have been proposed for the solution 
of problems of this type (see, e.g., [^3[ [3j, Q). The sim- 
plest procedure is the r approximation which was widely 
used for study of different problems of turbulent trans- 
port (see, e.g., |32|, |34|, [55|, [36|). One of the simplest pro- 
cedures that allows us to express the third moments /, 
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<f> N , . . . , H N in Eqs. flA15| ), flAl6| ) and flA23p in terms of 
the second moments, reads 



A N (k)-A^(k) = - 



A(k)-A'°)(k) 



(11) 



where the superscript (0) corresponds to the background 
turbulent convection (i.e., a turbulent convection with 
VjC/j = 0), and r(k) is the characteristic relaxation time 
of the statistical moments. Note that we applied the r- 
approximation ( [ll|) only to study the deviations from the 
background turbulent convection which are caused by the 
spatial derivatives of the mean velocity. The background 
turbulent convection is assumed to be known. 

The r-approximation is in general similar to Eddy 
Damped Quasi Normal Markowian (EDQNM) approx- 
imation. However there is a principle difference between 
these two approaches (see @ |§). The EDQNM clo- 
sures do not relax to the equilibrium, and this proce- 
dure does not describe properly the motions in the equi- 
librium state. Within the EDQNM theory, there is no 
dynamically determined relaxation time, and no slightly 
perturbed steady state can be approached j32] . In the r- 
approximation, the relaxation time for small departures 
from equilibrium is determined by the random motions in 
the equilibrium state, but not by the departure from equi- 
librium [p2[ . Analysis performed in |32j| showed that the 
r-approximation describes the relaxation to the equilib- 
rium state (the background turbulent convection) more 
accurately than the EDQNM approach. 

(c). We assumed that the characteristic times of vari- 



ation of the second moments (k) , (k) 



HQs.) are 



substantially larger than the correlation time r(k) for all 
turbulence scales. This allowed us to determine a station- 
ary solution for the second moments (k) , $^ (k) , . . . , 
H(k). 

(d). For the integration in k-space of the second mo- 
ments fij (k), $i (k), . . . , H(k) we have to specify a model 
for the background turbulent convection. Here we used 
the following model of the background turbulent convec- 
tion which is discussed in more details in Appendix B: 



/f(k) 
,(o) 



= /,[P w (k) + eP£->(kx)]W(*) 



*j 0) (k) = fc: 2 [fc 2 $(°)(k) ej P y -(k) 
+ip(°)(k)(exk) i ] , 



(12) 



(13) 
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d>(°)(k) = $: 



2a 



F<°>(k) = -6i/ (0) (k)(** • (exk)) 
G(°)(k) = (l + £ )/,/(°)(k)fc 2 , 
ff(°)(k) = 2H m W(k) , 



W(fe), (14) 

(15) 
(16) 
(17) 



e is the unit vec- 



where W(k) = W(k)/8i:k 2 , / (0) (k) = (fc ± /fe) 2 vj7(fc), 
e = (2/3)[(u^)/(u 2 ) — 2] is the degree of anisotropy of 
the turbulent velocity field u = u±+u z e, a is the degree 
of anisotropy of the turbulent flux of entropy (see below 
and Appendix B), Py(k) = <Jy — fey, fey = kikj/k 2 , 

k = k ± + k z e, k z = k • e, P£-\kj ) = 5 i4 - fe 

= (kx)i(kj_)j/fc^, ey = etej, 
tor directed along the z axis. Here r(fe) = 2rof(fe), 
M/(fc) = -df(k)/dk, f(fc) = {k/k o y-i, 1 < g < 3 is 
the exponent of the kinetic energy spectrum (q = 5/3 
for Kolmogorov spectrum), fco = and is the max- 
imum scale of turbulent motions, To = Iq/uq and uq is 
the characteristic turbulent velocity in the scale lo. Mo- 
tion in the background turbulent convection is assumed 
to be non-helical. In Eqs. ( p^ ) and jl^ ) we neglected 
small terms - 0(A/*;V/*) and ~ 0(A$*;V$*), re- 
spectively. Note that /; ( ° 3 (k)ey = /*/ (0) (k)- Now we 



calculate = J /^(k) dk using Eq. @ 



f (0) 



/ (0) = /. 



-Si 



^ (£y ey , 



(18) 



Note that = / $(°)(k)dk 
can be presented in the form 



The parameter a 



l+g(g + l)/(g-l) 

(Ix/l,) 9 ' 1 - 1 , 



(19) 
(20) 



where l± and Z z are the horizontal and vertical scales 
in which the correlation function $i° (r) = (s(x) u(x + 
r)) tends to zero (see Appendix B). The parameter £ 
describes the degree of thermal anisotropy. In particular, 
when l± = l z the parameter £ = and a = 1 . For ij_ <C / 2 
the parameter £ = — 1 and a — — 3/(g— 1). The maximum 
value £ max of the parameter £ is given by £ max = 9 — 1 for 
a = 3. Thus, for a < 1 the thermal structures have the 
form of column or thermal jets (l± < l z ), and for a > 1 
there exist the "pancake" thermal structures (l± > l z ) 
in the background turbulent convection. For statistically 
stationary small-scale turbulence the degree of anisotropy 
of turbulent velocity field varies in the range 



f 4(g + 3) 2(19 
1% 



a) 



1)' 
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< £ < OO 



(21) 



The negative (positive) degree of anisotropy e of a turbu- 
lent velocity field corresponds to that the vertical size of 
turbulent eddies in the background turbulent convection 
is larger (smaller) than the horizontal size. 



(e). In order to determine values /*, <&* and in the 
background turbulent convection we used balance equa- 
tions (|A5|)-( A7) for the second moments (see Appendix 
A). 



C. Turbulent flux of entropy 

The procedure described in this Section allows us to 
determine the Reynolds stresses and turb ulent flux of 
entropy which are given by Eqs. ( |A33| ) and (A34) in Ap- 
pendix A, where we considered the case V ■ U = 0. 
In particular, the formula for turbulent flux of entropy 
reads: 



-5a(V • U±)*; + {a + 3/2)(u>x$;) 



+3(u> || x<fr*)] 
+ (E- V)U, 
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(VxT) 



(22) 



where 



(VxU) z , *; = $*e, 

(2/5)r (<7 - 2)(** • (exU)) and E 
2a(g + 1)/3]<I>* - 3**} 



(l/5)r {[2 - q - 
It is shown below that the first 



and the second terms in Eq. (22) are responsible for the 
large-scale convective wind instability in a shear-free tur- 
bulent convection (see Section IV), while the third term 
in the turbulent flux of entropy (]22|) causes the convective 
shear instability (see Section V). 

The turbulent flux of entropy can be obtained even 
from simple symmetry reasoning. Indeed, this flux can be 
presented as a sum of two terms: (su) = $* + Pijk^jUk, 
where <&* determines the contribution of the Kolmogorov 
turbulence and it is independent of ViUj, whereas the 
second term is proportional to WiUj and describes the 
contribution of the tangling turbulence. Here 0ijk is an 
arbitrary true tensor and U is the mean velocity. Using 
the identity Vjt/j = (<5U)y — {l/2)eijkOJk, the turbulent 
flux of entropy becomes 

(sm) = <f>* + rnjQj + {u>x5)i + //,,/. (rill) , (23) 



where (6XJ)i 



(ViUj + VjUi)/2, uj = VxU is the 



mean vorticity, and £yfc is the fully antisymmetric Levi 
Civita tensor. In Eq. (23), rjij is a symmetric pseudo- 
tensor, 8 is a true vector, fi^jk is a true tensor symmet 
ric in the last two indexes, $ 
vectors. The tensors rjij, fiijk 
constructed using two vectors: 
vector e. For example, rjij — 
and \iijk = A 3 $* z e ijk + j4 4 $ 
known coefficients and ey^ = 
lowing expression for the turbulent flux of entropy in a 
divergence-free mean velocity field: 

* = **-(a 3 + a 4 )(v-u ± )*; 

+ {A X + A 2 )(u)x<f>*) + Ai(uj {1 x**) 
-A±{V ■ Ul)#I] , (24) 



= (su) and <&* are true 
and the vector 5 can be 
<&* and the vertical unit 
0, 6 = A^* + A 2 $*e, 
■ejk, where Ak are the un- 
: eiCjek- This yields the fol- 
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where U = Uj_ + U z e, 3?* = $*e and u> y = u z e. Equa- 
tions ( p2[ ) and coincide if one sets A\ = To(q + 1)/5, 
A 2 = toI<7 + l)(a - 3/2)/15, A3 = T a(q + l)/3, and 
At = 0. Note that ** = -k t VS - T ^*(dU^>(z)/dz), 
U(°) (z) is the imposed horizontal large-scale flow velocity 
(e.g., a wind velocity). 



III. CONVECTIVE WIND INSTABILITY IN A 
SHEAR-FREE TURBULENT CONVECTION 

In this section we studied the mean-field dynam- 
ics for a shear-free turbulent convection. We showed 
that under certain conditions a large-scale instability 
is excited, which causes formation of large-scale semi- 
organized structures in a turbulent convection. 

The mean-field dynamics is determined by Eqs. (|]) and 
(||) . To study the linear stage of an instability we derived 
linearized equations for the small perturbations from the 
equilibrium^ UP = U z - U { z cq \ u>W = q _ ^(eq) an d 
5(1) = S - : 



A 



dU. 



(i) 



at 
at 



d_ 

dz 



-<7A±S (i: 
■(exV)i 



v,/ (1) 

v 3 J 13 



dU, 



(eq) 



0t 



iV/, 



dz 



(25) 

r/W ,(26) 
{/« ,(27) 



where /, 



(i) 



fij L 



(0) 



and the Reynolds stresses /y are 
given by Eqs. ( A33 ) in Appendix A, Aj_ = A — d 2 /dz 2 
N b = N b e and 

V-*W = _( To /30)( g + 1){(** ■e)[10aA ± 

-(8a-3)ApP +6((**xe) • V)u)W} 

-KyViVi^ , (28) 

«« = r ,544 0) + (l/2) 5 To(5,(4-7)(e,:$* 

+ej*i)] • (29) 

Equation (E~~ 



follows from Eqs. (A6) and (A7) 



A. The growth rate of convective wind instability 



Let us consider a shear-free turbulent convection 



(ViJ7. 

S (eq) 



(0) 



0) with a given vertical flux of entropy 

§ y z q> e. We also consider an isentropic basic reference 
state, i.e., we neglect terms which are proportional to 
(N b + dS (c ^ /dz)U^ ] in Eq. @. We seek for a solution 
of Eqs. (gg)-(gg) in the form oc exp(7i nst t+iK-R), where 
K is the wave vector of small perturbations and 7i ns t is 
the growth rate of the instability. Thus, the growth rate 
of the instability is given by 



7inst 



v T K 2 Ay\ + AB/A 2 - l]/2 , 



(30) 



where 

A = B 1 + B 2 , B = f3X{c 7 -c 8 X)-B 1 B 2 , (31) 

B 1 = c l + c 6 X - c 3 X 2 , B 2 = c 4 - c 5 A, ci = {q + 3)/5, 
c 3 = e(q + 3)/4, c 4 = S„(2 + 3a), c 5 = 36, (a - e/2), 
cq = e(q + 5)/4, C7 = //(8a — 3)/10, eg = ya, with 
a = a*(4 - 7 )(1 + e/2), y, = 6a*(g + 1)(1 + e/2)/5„ 
= (I0K)- 2 , X = sin 2 9, a* = 2S^ ( z q) gr /f*, and 6 
is the angle between e and the wave vector K of small 
perturbations. Here we used that in equilibrium $i cq ^ = 
When ^ > 1 the growth rate of the instability is 
given by 



7inst oc v T K 2 v^lsinf 



3 5a 2 
a--- T sm 



-1 1/2 



(32) 



Thus for large the growth rate of the instability is pro- 
portional to the wave number K and the instability oc- 
curs when a(5cos 2 6 — 1) > 3/2. This yields two ranges 
for the instability: 



2(5 cos 2 6 - 1) 
3 



< a < 3 , 



9-1 



< a < - 



2(1 - 5 cos 2 1 



(33) 
(34) 



where we took into account that the parameter a varies 
in the interval —3/{q — 1) < a < 3 (see Appendix B). 
The first range for the instability in Eq. ( |33"| ) is for the 
angles 3/10 < cos 2 9 < 1 (for q — 5/3, the aspect ra- 
tio < L z /L± < 1.53), and the second range (|34l) for 
the instability corresponds to the angles < co?6* < 
(3 — q)/10 (the aspect ratio 2.55 < L z /L± < 00), where 
L z /L± = K ± /K z = tan 6. The conditions (§§) and (§4|) 
correspond to V • <&W < 0. 

Figure [l] demonstrates the range of parameters L z /L± 
and L/lo where the instability is excited, for different 
values of the parameter a (from —4.5 to 3) and dif- 
ferent values of the parameter e = —1; 0; 5. Here 

L = 1/ V/ ' L z 2 + LJ 2 and we assumed that a* = 1. The 
threshold of the instability L cr depends on the parameter 
e. For example, for a — 3 the threshold of the instabil- 
ity L cl - varies from 31q to 71q (when s changes from — 1 
to 5). The negative (positive) degree of anisotropy e of 
turbulent velocity field corresponds to that the vertical 
size of turbulent eddies in the background turbulent con- 
vection is larger (smaller) than the horizontal size. The 
reason for the increase of the range of instability with 
the decrease of the degree of anisotropy e is that the rate 
of dissipation of the kinetic energy of the mean velocity 
field decreases with decrease of e and it causes decrease 
of the threshold of the instability. The instability does 
not occur when 1.53 < L z /L± < 2.55 for all e. 

Figure || shows the growth rate of the instability as 
function of the parameter L/Iq (FIG. 2a) and of the pa- 
rameter L z /L ± (FIG. 2b) for e = and a = 2 (the 
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FIG. 1: The range of parameters (L z /L±; L/Iq) for which the 
convective wind instability occurs, for different values of the 
parameter a: (from —4.5 to 3) and for different values of the 
parameter e: a), e = 5 ; b). e — ; c). e = — 1. 



FIG. 2: The growth rate of the convective wind instability 
as functions of: a). L/lo (for different values of parameter 
L z /L± = 0.3; 1; 1.23); and b). L 2 /L± (for different values of 
parameter L/lo = 4.2; 4.6; 9.6; 20) for e = and a = 2. 



first range of the instability). This range of the insta- 
bility corresponds to the "pancake" thermal structures 
of the background turbulent convection (l z /l± « 2/3 for 
a = 2). The maximum of the growth rate of the insta- 
bility (7 m ax ~ 0.045 Tq 1 ) reaches at the scale of pertur- 
bations L m w 9.4 ?o (for L z /L± « 0.76). In this case the 
threshold of the instability L cr 4.2 Zo- 

Figure [5] demonstrates the growth rate for the second 
range of the instability (a = —3). Note that this range 
of the instability corresponds to the thermal structures 
of the background turbulent convection in the form of 
columns (l z /l± » 2 for a = — 3). In contrast to the first 
range of the instability, the growth rate increases with 
L z /L± in the whole second range of the instability (see 
FIG. 3b). 



B. Mechanism of the convective wind instability 

The convective wind instability results in formation 
of large-scale semi-organized structures in the form of 
cells (convective wind) in turbulent convection. The 
mechanism of the convective wind instability, associated 
with the first term €> oc — roo;(V • U^)$* in the ex- 
pression for the turbulent flux of entropy [see Eq. (p2|)], 
in the shear-free turbulent convection at a > is as 
follows. Perturbations of the vertical velocity U z with 
dU z /dz > have negative divergence of the horizontal 
velocity , i.e., divU^ < (provided that divU ps 0). 
This results in the vertical turbulent flux of entropy 
$2 oc — divUj_, and it causes an increase of the mean en- 
tropy (S cx p- 1 ' 2 U z $*/ul) [see Eqs. and ©]. 

On the other hand, the increase of the the mean en- 
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FIG. 3: The growth rate of the convective wind instability 
as functions of: a). L/lo (for different values of parameter 
L z /L± = 3.5; 4; 5); and b). L z /L± (for different values of 
parameter L/lo = 6; 10; 20; 40); for e = 5 and a = —3. 



tropy increases the buoyancy force oc gS and results 
in the increase of the vertical velocity U z oc tq/3 1 / 2 gS 
and excitation of the large-scale instability [see Eqs. ( |2q ) 
and (|32|)]. Similar phenomenon occurs in the regions 
with dll z /dz < whereby divUj^ > 0. This causes a 
downward flux of the entropy and the decrease of the 
mean entropy. The latter enhances the downward flow 
and results in the instability which causes formation of 
a large-scale semi-organized convective wind structure. 
Thus, nonzero divU^ causes redistribution of the ver- 
tical turbulent flux of entropy and formation of regions 
with large vertical fluxes of entropy (see FIG. [|). This 
results in a formation of a large-scale circulation of the 
velocity field. This mechanism determines the first range 
for the instability. 

The large-scale circulation of the velocity field causes a 
nonzero mean vorticity Q, and the second term [propor- 




FIG. 4: The effect of a nonzero div which causes a re- 
distribution of the vertical turbulent flux of the entropy and 
results in a formation of a large-scale circulation of the ve- 
locity field. Fluid flow with divtjj_ < (a) produces regions 
with vertical fluxes of entropy and vertical fluid flow (b) in 
these regions. 



tional to (a + 3/2)(wx$*)] in the turbulent flux of the 
entropy (|2^) is responsible for a formation of a horizon- 
tal turbulent flux of the entropy. This causes a decrease 
of the growth rate of the convective wind instability (for 
a > 0), because it decreases the mean entropy S in the 
regions with dU z /dz > 0. The net effect is determined by 
a competition between these effects which are described 
by the first and the second terms in the turbulent flux of 
the entropy ((H) . The latter determines a lower positive 
limit a m in = 3/8 of the parameter a. 

When a < —3/2 the signs of the first and second 
terms in the expression ( p2| ) for the turbulent flux of en- 
tropy change. Thus, another mechanism of the convec- 
tive wind instability is associated with the second term 
in the expression ( p2| ) for the turbulent flux of entropy 
when a < —3/2. This term describes the horizontal flux 
of the mean entropy <& y oc r (a + 3/2)(u>X<&*). The lat- 
ter results in the increase (decrease) of the mean entropy 
in the regions with upward (downward) fluid flows (see 



8 




FIG. 5: The effect of a nonzero wxtj which induces the hor- 
izontal flux of the mean entropy & y and causes increase (de- 
crease) the mean entropy in the regions with upward (down- 
ward) fluid flow when a < —3/2. 

FIG. ||). On the other hand, the increase of the mean 
entropy results in the increase of the buoyancy force, the 
mean vertical velocity U z and the mean vorticity ui. The 
latter amplifies the horizontal turbulent flux of entropy 
3> y and causes the large-scale convective wind instability. 
This mechanism determines the second range for the con- 
vective wind instability. The first term in the turbulent 
flux of entropy at a < causes a decrease of the growth 
rate of the instability because, when a < 0, it implies a 
downward turbulent flux of entropy in the upward flow. 
This decreases both, the mean entropy and the buoyancy 
force. Note that, when a < —3/2, the thermal structure 
of the background turbulence has the form of a thermal 
column or jets: l z /l± > 3.34. Even for a < 0, the ratio 
l m /l x > 1-54. 



IV. CONVECTIVE SHEAR INSTABILITY 

Let us consider turbulent convection with a linear shear 
tj(°J = (X/To)ze y and a nonzero vertical flux of en- 
tropy = <&z ^ e, where A is dimensionless parame- 
ter which characterizes the shear. We also consider an 
isentropic basic reference state, i.e., we neglected a term 
which is proportional to {N b +dS {c ^ Jdz)ui 1] in Eq. @. 
We seek for a solution of Eqs. (p5|)-(p7|) in the form 
tjW = V exp(7i nst t) cos(f2i + K • R). Here, for simplicity, 
we study the case K • Ij( c< i) = 0. 



A. The growth rate of convective shear instability 

Using a procedure similar to that employed for the 
analysis of the convective wind instability we found that 
the growth rate of the convective shear instability is de- 
termined by a cubic equation 



(7 + S 3 )(7 2 + Aj - B) + 8/3 2 7o 3 = , 



(35) 



where 7 = (7 inst 

Cg = 18/i6*/5, i 
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FIG. 6: The range of parameters (L z /L±; L/Iq) for which the 
convective shear instability occurs, for a = 2, e = and a = 
—3, e — 5 for different values of the shear A = 0.05; 0.1; 0.15. 



B3 = c\ + C2X, C2 = e{q + l)/4. The growth rate of 
the instability for (3 3> 1 is given by 



7inst ~ v T K z ( (3 2 ' z 70 + P 



3l/3 



c 

1270" 



(36) 



where C — X{c-j — c%X). The instability results in gen- 
eration of the convective shear waves with the frequency 



n ~ VSvtK 2 ( (3 2/3 70 - /3 1/3 



C 
1270" 



(37) 



The flow in the convective shear wave has a nonzero hy- 
drodynamic helicity 



X = V- (VxV) = - 



2A^V 2 

ro(n 2 + 7Lt) 



(38) 



-4 eq) (l + e/2)/(A$i cq) ) and 



Therefore, for A > the mode with K x > has a negative 
helicity and the mode with K x < has a positive helicity. 

Figure ^| shows the range of parameters L z /Lj_ and 
L/Iq where the convective shear instability occurs, for 
a = 2, e = and for different values of the shear A = 
0.05; 0.1; 0.15. There are two ranges for the instability. 
However, even a small shear causes an overlapping of the 
two ranges for the instability and the increase of shear 
(A) promotes the convective shear instability. 

Figures fj] and || demonstrate the growth rates of the 
convective shear instability and the frequencies of the 
generated convective shear waves for the first (a — 2) and 
the second (a = —3) ranges of the instability. The curves 
in FIGS. |^-|| have a point L* whereby the first derivative 
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d"finst/dK has a singularity. At this point there is a bifur- 
cation which is illustrated in FIGS. 7] and |[ The growth 
rate of the convective shear instability is determined by 
the cubic algebraic equation ( |35| ) . Before the bifurcation 
point [L < L„) the cubic equation has three real roots 
(which corresponds to aperiodic instability). After the 
bifurcation point (L > the cubic equation has one 
real and two complex conjugate roots. When L > L* 
the convective shear waves are generated. When the pa- 
rameter L z /L± increases, the value decreases. When 
L z > L±, the bifurcation point L» < L CI . For a given pa- 
rameter L/Iq there are the lower and the upper bounds 
for the parameter L z /L± when the convective shear in- 
stability occurs. For large enough parameter L — L u the 
upper limit of the range of the instability does not ex- 
ist, e.g., for A = 0.05 the parameter L u = 47?o and for 
A = 0.15 the parameter L u = 13/o- 

Note that when L < L» the convective shear waves are 
not generated and the properties of the convective shear 
instability are similar to that of the convective wind in- 
stability (compare FIG. 2b and the curve for L/Iq = 6 
in FIG. 8c). However for L > L* these two instabilities 
are totally different. The properties of the convective 
shear instability in the first and in the second ranges of 
the instability are different. In particular, in the second 
range of the convective shear instability the growth rate 
monotonically increases, and the frequency of the gener- 
ated convective shear waves decreases with the parameter 
L,/L ± . 



B. Mechanism of convective shear instability 

The mechanism of the convective shear instabi lity 
associated with the last term in the expression ( |22] ) 
for the turbulent flux of entropy [<fr oc na(u)„x&*)] 
is as follows. The vorticity perturbations u> = 
(VxU) z generate perturbations of entropy: S oc 
/3~ 1//6 (tq/u )$*\~ 2 / 3 exp(i7r/6)u). Indeed, consider two 
vortices (say, "a" and "b" in FIG. ||) with the opposite 
directions of the vorticity u> » . The turbulent flux of en- 
tropy is directed towards the boundary between the vor- 
tices. The latter increases the mean entropy between the 
vortices ("a" and "b"). 

Similarly the mean entropy between the vorticies "b" 
and "c" decreases (see FIG. ||). Such redistribution 
of the mean entropy causes increase (decrease) of the 
buoyancy force and formation of upward (downward) 
flows between the vortices "a" and "b" ("b" and "c"): 
Ug cx /3 1 / 3 A _2/ ' 3 5ro exp(— m/2>) S. Finally, the vertical 
flows generate vorticity Co cx — /3 -1 / 6 A 1 / 3 exp(i7r/6)J7 z /^o, 
etc. This results in the excitation of the instability with 
the growth rate 7; ns t oc K 2 ! 3 and generation of the con- 
vective shear waves with the frequency fl oc K 2 I 3 . For 
perturbations with K x — the convective shear insta- 
bility does nor occur. However, for these perturbations 
with K x = the convective wind instability can be ex- 
cited (see Section III), and it is not accompanied by the 
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FIG. 7: The growth rates of the convective shear instability 
[a) and c)] and the frequencies of the generated convective 
shear waves [b) and d)] for the first (a = 2) range of the 
instability and for e = 0. Corresponding dependencies on the 
parameter L/Iq are given for different L z /L± and visa versa. 



generation of the convective shear waves. We considered 
here a linear shear for simplicity. The equilibrium is also 
possible for a quadratic shear, i.e., when tj(°) = Xz 2 e y . 



V. DISCUSSION 

The "convective wind theory" of turbulent sheared 
convection is proposed. The developed theory predicts 
the convective wind instability in a shear-free turbulent 
convection. This instability causes formation of large- 
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scale semi-organized fluid motions (convective wind) in 
the form of cells. Spatial characteristics of these motions, 
such as the minimum size of the growing perturbations 
and the size of perturbations with the maximum growth 
rate, are determined. 

This study predicts also the existence of the convec- 
tive shear instability in the sheared turbulent convec- 
tion. This instability causes formation of large-scale 
semi-organized fluid motions in the form of rolls (some- 
times visualized as the boundary layer cloud streets). 



FIG. 9: The effect of a nonzero On x<]?* which causes a re- 
distribution of the horizontal turbulent flux of the entropy. 
For two vortices ("a" and "b") with opposite directions of the 
vorticity £5.. , the turbulent flux of entropy is directed towards 
the boundary between the vortices. The latter increases the 
mean entropy between the vortices ("a" and "b"). Similarly, 
the mean entropy between the vorticies "b" and "c" decreases. 



These motions can exist in the form of generated con- 
vective shear waves, which have a nonzero hydrodynamic 
helicity. Increase of shear promotes excitation of the con- 
vective shear instability. 

The proposed here theory of turbulent sheared convec- 
tion distinguishes between the "true turbulence", corre- 
sponding to the small-scale part of the spectrum, and 
the "convective wind" comprising of large-scale semi- 
organized motions caused by the inverse energy cas- 
cade through large-scale instabilities. The true turbu- 
lence in its turn consists of the two parts: the famil- 
iar " Kolmogorov-cascade turbulence" and an essentially 
anisotropic "tangling turbulence" caused by tangling of 
the mean-velocity gradients with the Kolmogorov-type 
turbulence. These two types of turbulent motions over- 
lap in the maximum-scale part of the spectrum. The 
tangling turbulence does not exhibit any direct energy 
cascade. 

It was demonstrated here that the characteristic length 
and time scales of the convective wind motions are much 
larger than the true-turbulence scales. This justifies sep- 
aration of scales which is required for the existence of 
these two types of motions. It is proposed that the term 
turbulence (or true turbulence) be kept only for the Kol- 
mogorov and tangling turbulence part of the spectrum. 
This concept implies that the convective wind (as well 
as semi-organized motions in other very high Reynolds 
number flows) should not be confused with the true tur- 
bulence. The diagram of interactions between turbulent 
and mean flow objects which cause the large-scale in- 
stability and formation of semi-organized structures is 
shown in FIG. |h]. 

Now let us compare the obtained results with the prop- 
erties of semi-organized structures observed in the atmo- 
spheric convective boundary layer. The semi-organized 
structures are observed in the form of rolls (cloud streets) 
or three-dimensional convective cells (cloud cells). Rolls 
usually align along or at angles of up to 10° with the mean 
horizontal wind of the convective layer, with lengths from 
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flow objects which cause a large-scale instability. 



20 to 200 km, widths from 2 to 10 km, and convective 
depths from 2 to 3 km |2| . The typical value of the aspect 
ratio L z /L± w 0.14— 1. The ratio of the minimal size of 
structures to the maximum scale of turbulent motions 
is L/Iq = 10 — 100. The characteristic life time of rolls 
varies from 1 to 72 hours [0 . Rolls may occur over both, 
water surface and land surfaces. The suggested theory 
predicts the following parameters of the convective rolls: 
the aspect ratio L z /Lj_ ranges from very small to 1, and 
L/Iq = 10 — 100. The characteristic time of formation of 
the rolls ~ ro/7i ns t varies from 1 to 3 hours. The life 
time of the convective rolls is determined by a nonlinear 
evolution of the convective shear instability. The latter 
is a subject of a separate ongoing study. 

Convective cells may be divided into two types: open 
and closed. Open-cell circulation has downward motion 
and clear sky in the cell center, surrounded by cloud 
associated with upward motion. Closed cells have the 
opposite circulation Both types of cells have di- 

ameters ranging from 10 to 40 km and aspect ratios 
L z /Lj_ « 0.05— 1, and both occur in a convective bound- 
ary layer with a depth of about 1 to 3 km. The ra- 
tio of the minimum size of structures to the maximum 
scale of turbulent motions is L/Iq = 5 — 20. The de- 
veloped theory predicts the following parameters of the 
convective cells: the aspect ratio L z / L± ranges from very 
small to 1, and L/Iq = 5 — 15. The characteristic time of 



formation of the convective cells ~ To/7i ns t varies from 
1 to 3 hours. Therefore the predictions of the devel- 
oped theory are in a good agreement with observations 
of the semi-organized structures in the atmospheric con- 
vective boundary layer. Moreover, the typical temporal 
and spatial scales of structures are always much larger 
then the turbulence scales. This justifies separation of 
scales which was assumed in the suggested in the theory. 
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APPENDIX A: DERIVATIONS OF 
EXPRESSIONS FOR THE REYNOLDS STRESSES 
AND TURBULENT FLUX OF ENTROPY 

Equations (||) and (|J) yield the following conservation 
equations for the kinetic energy W v = poV 2 /2, for Ws = 
p S 2 /2 and for W* = p Sv : 



dW v 

at 

dW s 

~dT ' 
at 



■ V-F„ 
V-F s 



V,F? 



Is-D s , 



Df 



(Al) 
(A2) 
(A3) 



where I v = -p (v • g)S, I s = -I v flf/g 2 , and I* = 
-p [v(v • N fc ) + 5 2 g] + (P/p )V(S>o) are the source 
terms in these equations, D v = — p (v ■ f„), Dg = 
p S(V • F„) and D* = -p Sf v + (p /T )v(V • F„) are 
the dissipative terms, F„ = v(W„ + P), F$ — vWs 
and F* = poS vjVj + SPSij are the fluxes. Equations 



( |Al| ) and (A2) yield conservation equation for We = 
WM 2 b /g 2 + W S 



dW E 



V-F, 



-D T 



(A4) 



where De = D v ^l 2 /g 2 + Ds is the dissipat ive term, and 
F E = Fyflf/g 2 + F s is the flux. Equation (|A^ ) does not 



have a source term, and this implies that without dissi- 
pation (De — 0) the value J We dV is conserved, where 
in the latter formula the integration is performed over 
the volume. For the convection Q 2 < and, therefore, 

w s ~w v \nl\/g 2 . 
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Using Eqs. (Al)-( A3) we derived balance equatio ns fo r 
the second moments. In particular, averaging Eqs. ( |Al|) - 
(A3) over the ensemble of fluctuations and subtracting 
from these equations the corresponding equations for the 
mean fields: poU 2 /2, p S 2 /2, PqSIJ, yields 



d_ 

dt 



0_ 

dt 



U V / pp + 2/yVjUi + 2*-. 



-V 



Po 



J*_ 

Tq<5* 



(A5) 



U • V + • V)U t + fij (N b + VS%- 

Po " ' \ 



~ gei (4-7)H+— V^ij- [ — )/i(N b !, 
2 Po \ Po / 



d_ 

dt 



2ro<5* ^ ^ ' 
U V )P + 2$ i (N b + V5)j 



tq<5* 



(A6) 
. (A7) 



where Pr = ^/k is the Prandtl number, v is the kinematic 
viscosity, * B = p f pp u/2 + (pu), *y = p {suiUj) + 
5ij('yPoH/2 + (sp)), and we took into account that the 
dissipations of energy, the flux of entropy and the second 
moment of entropy H are determined by the background 
turbulent convection described by Eqs. (|T^)-([l7j). In 
derivation of Eq. ( |A6| ) we used an identity pqSV(P/ po) — 
V( 7 P S 2 )/2 + 7^oS 2 N b + (7/2 - l)p gS 2 , and we as- 
sumed that S ~ P/7P0, » .e. ; we ne glected fluctuations of 
density p/po- Equations (|A5|)-(A7) allow us to determine 
/*, and -ff* in the background turbulent convection 
(see below). 

Using Eqs. (j(|) and ([?]) we derived equations for the 
following second moments: 



Ai(k,R) = 


(Mi(k,R)^(-k,R)) , 


(A8) 


$ 4 (k,R) = 


(s(k,RK(-k,R)) , 


(A9) 


F(k,R) = 


(s(k,R)w(-k,R)) , 


(A10) 


G(k,R) = 


(w(k,R)w(-k,R)) , 


(All) 


ff(k,R) = 


( s (k,R)s(-k,R)) , 


(A12) 



where oj = (Vxu) z and we use a two-scale approach, , 
i.e., a correlation function is written as follows 

(ui(x)uj(y)) = J {ui(ki)uj(k.2)) exp[i(k r x 
+k 2 -y)]cik 1 dk 2 



/ii(k,R) = 



J fij (k, R) exp(ik-r) dk , 
/^(k + K/2)^(-k + K/2)) 



x exp(iK-R) dK 



(see, e.g., 37, f38|), where R and K correspond to the 
large scales, and r and k to the small scales, i.e., R = 
(x + y)/2, r = x-y, K = ki + k 2 , k = (ki - 



k 2 )/2. This implies that we assumed that there exists a 
separation of scales, i.e., the maximum scale of turbulent 
motions Iq is much smaller then the characteristic scale L 
of inhomogeneities of the mean fields. In particular, this 
implies that r < Iq <C R. Our final results showed that 
this assumption is indeed valid. Now let us calculate 



df l] (k 1 ,k 2 ) 
dt 



a$j(ki,k 2 ) 



du n (k 



dt 



-u 3 (k 2 )) 



+ K(k 1 )P jn (k 2 )^ 2 l),(A13) 
(s{k 1 )P jn (k 2 ) — ) 



at 



v at 



Uj(k 2 )) 



(A14) 



where we multiplied equation of motion (^) rewritten in 
k-space by Pjj (k) = <5jj — fcy in order to exclude the 
pressure term from the equation of motion. 

Thus, equations for /y(k, R) and 3?(k, R) read: 



a/«j(k) 
at 

dt 



— Iijmnfmn(k) ~\- Nij(k) , (A15) 

= 4^(k)+Mi(k) , (A16) 



where 



L ijmn 



J iq w jn u mp 



^im^jq^np 



_d_ 

l dK 



(A17) 



5 



M<(k) 



2^n^ jU n -\- $ijk n 

OKm 

ge m [P lm (k)^(k) + P jm (k)$ i (-k)] 
+/#(k), 



VmUn-VjUi , (A18) 

(A19) 



■/™(Nb + V5) f 
-4>f , 



ge m Pim(k)H 



(A20) 



and hereafter we consider the case with V • U = (i.e., 
A = 0). Here /y and $f are the third moment s ap - 
peari ng d ue to the nonlinear terms. Equations ( |A15| ) 
and (A16) are written in a frame m ovin g with a local 
velocity U of the mean flow. In Eqs. ([A15|)-(|A20|) 



glected small terms which are of t he or d er of Q (V 3 U) and 
0(V 2 fif, V 2 $ t ). Note that Eqs. (|Al5j)-(|A20|) do not con- 
tain the terms prop ortio nal to Q(V 2 U). The first term 
in the RHS of Eqs. ( A15 ) and ( A16 ) depends on the gra- 
dients of the mean fluid velocity (Vjjj). Equations for 
the second moments G(k), P(k) and P(k) read: 



dG(k) - - . df l3 (k) 

= (exki)i(exk 2 )j 



dt 
dF(k) 

dH(k) 
~dt~ 



dt 



—i(exki)j 

Q(k) , 



dt 



(A21) 
(A22) 
(A23) 
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where Q(k) = -2*(k) • (N b + VS) + H N , and H N is the 
third moment appearing due to the nonlinear terms, ki — 
k-(i/2)V, k 2 = k+(z /2) V. The terms ~ * in the tensor 
A^-(k) [see Eqs. jAiq ) and ( |A"l9| )] can be considered as a 
stirring force for the turbulent convection. On the other 
hand, the terms ~ (N b + VS) in Eqs. (|A16fl, QA20D and 



( |A23j ) are the sources of the flux of entropy 3? and the 
second moment of entropy H. Note that a stirring force 
in the Navier-Stokes turbulence is an external parameter. 

Since the equations for the second moments contain 
the third moments, a problem of closure for the higher 
moments arises. In this study we used the r approx- 
imation [see Eq. (fl|)] which allows us to exp ress the 



third moments f£j, and H N in Eqs. (|A15[ ), ( |A16| ) 



and (A23) in terms of the second moments. Here we de- 
fine a background turbulent convection as the turbulent 
convection with zero gradients of the mean fluid veloc- 
ity (VjC/j = 0)]. The background turbulent convection is 
determined by the following equations: 



Ot 



M^(k) 



Ot 



QW(k) 



(A24) 
(A25) 
(A26) 



A nonzero gradient of the mean fluid velocity results 
in deviations from the background turbulent convection. 
These deviations are determined by the following equa- 
tions: 



fij fi 



Ot 

Ot T(fc) . 

d(H-H^) _ H-H<® 



(0) 



r(fc) 

(0) 



at 



/ m n(k) - "_„,; ,J ,(A27) 
(A28) 
(A29) 



r(fc) ' 



where the deviations (caused by a nonzero gradients of 
the mean fluid velocity) of the functions N t] (k) - iV^ 0) (k) 
and Mi(k) — (k) from the background state are de- 
scribed by the relaxation terms: —{fij — /A- )/r{k) and 
— ($j — $|°' ) )/T(fc), respectively. Similarly, the devia- 
tion Q(k) - Q (0) (k) is described by the term — {H — 
H^)/r{k). Here we assumed that the correlation time 
r{k) is independent of the gradients of the mean fluid 
velocity. 

Now we assume that the characteristic times of varia- 
tion of the second moments /y(k), $i(k) and H(k) are 
substantially larger than the correlation time r{k) for all 
turbulence scales. This allows us to determine a station- 
ary solution for the second moments /y(k), $i(k) and 
H(k) : 



/y(k) = /^(k)+r(fe)J, 



(A30) 



*i(k) 
J?(k) 



*i 0) (k)+r(fc)4-$f(k 
ff(°)(k), 



(A31) 
(A32) 



where we neglected the third and higher order spatial 
derivatives of the mean velocity field U. 

For the integration in k-space of the second moments 
fijQc), $i(k), . . . , H(k, R) we have to specify a model for 
the background turbulent convection. We used the model 
of the background turbulent convection determined by 
Eqs. (p^)-(p^|). For the integration in k-space we used 
identities give n in Appe ndix C. The integration in k- 
space of Eqs. ( |A30 ) and ( A31 ) yields the following equa- 
tions for the Reynolds stresses and the turbulent flux of 
entropy: 



fi. 



-ejVi)U z 





dz 



VjUi) + a 2 (eiVj 
ejUi) 



+ -^(c 3 e tj + a 3 dij) 



(A33) 



TO 

' 5 



2(**xe)w + 5(** • V)U 
(3-j)(exV)(** •(exU)) 
(?-l)[((**xe).V)(exU) 



b 2 U)- 
- 2($ : 



b 3 V^U Z 



V,)U 



(A34) 



where v T = t /*/6, uj = (VxU) 2 , a\ = c\ + c 2 , a 2 = 
-e(« - l)/4, a 3 = -e(5 - <?)/4, 6 X = (8a - 3)(g + 1) 
62 = 3(9 -q)- 2a(q + 1), b 3 = {2a + 3){q + 1), d = 
(q + 3)/5, c 2 = efa + D/4, eg - e(« + 3)/4. 



Equations ( A33 ) and (A34) imply that there are two 
contributions to the Reynolds stresses and turbulent flux 
of entropy which correspond to two kinds of fluctuations 
of the velocity field. The first contribution is due to 
the Kolmogorov turbulence with the spectrum (oc k~ q ), 
and it corresponds to the background turbulent convec- 
tion. The second kind of fluctuations depends on gradi- 
ents of the mean velocity field and is caused by a "tan- 
gling" of gradients of the mean velocity field by turbu- 
lent motions. The spectrum of the tangli ng tu rbulence 
is W{k)r{k) cx fc 1 " 29 [see Eqs. ( |A30| ) and (|A3l|)]. These 
fluctuations describe deviations from the background tur- 
bulent convection caused by the gradients of the mean 
fluid velocity field. 

Now we calculate a dissipation of the kinetic energy of 
the mean flow U : 



D V = -(l/2)(/« - fffWiUj + VjUt) 



(A35) 



using a general form of the velocity field f/j = 
Vi{t, K) exp(iK ■ R), where 

Vi{t,K) = f-Q [P l3 {K)e 3 V z {t,K) 

-iK- 2 (exK)iu(t,K)], (A36) 
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and Q = (VxV) z . The result is given by 



D,U. K) }>r \\ n [jL) [lif 2 i?(t,K)+£ 2 (/.K) 



-6 5 if 2 K 2 (i,K) , 



(A37) 



where 64 = Ci + c 2 sin 2 and b 5 = a 2 + c 3 cos 2 9. The 
function Du(t,K) must be positive for statistically sta- 
tionary small-scale turbulence. The latter is valid when 
e satisfies condition (|2l|). 



Here we also used the identity (kixej^kixejj 



k 2 P 
tain 



(->-)/ 



(k±). Substituting Eq. (|B3j) into Eq. (jBJ) we ob 



(fc/fex) 2 {/*/ (0) (k)i^-(k) 
-/,/(°)(k)]^|(k ± )}. 



[G^(k)/fc 2 



(B4) 



Thus two independent functions determine the corre- 
lation function of the anisotropic turbulent velocity 
field. In the isotropic three-dimensional turbulent flow 
G'°'(k)/fc 2 = /»/(°)(k) and the correlation function 
reads 



APPENDIX B: THE MODEL OF THE 
BACKGROUND TURBULENT CONVECTION 

A simple approximate model for the three-dimensional 
isotropic Navier-Stokes turbulence is described by a 
two-point correlation function of the velocity field 
fij(t,x,y) = {ui(t,x)uj(t,y)) with the Kolmogorov spec- 
trum W(k) oc k~ q and q = 5/3. The turbulent con- 
vection is determined not only by the turbulent veloc- 
ity field u(t, x) but by the fluctuations of the entropy 
s(i,x). This implies that for the description of the tur- 
bulent convection one needs additional correlation func- 
tions, e.g., the turbulent flux of entropy <&j(t, x, y) = 
(s(t, x)Ui(t, y)) and the second moment of the entropy 
fluctuations H(t,x.,y) — (s(t, x)s(t, y)). Note also that 
the turbulent convection is anisotropic. 

Let us derive Eqs. ([[2]) and (|l3|) for the correlation 
functions and To this end, the velocity is 
written as a sum of the vortical and the potential compo- 
nents, i.e., u± — Vx(Ge) + Vj_0, where ui = (Vxu) z = 
-A X C, A ± = -du z /dz, V_l = V-e(e- V). Hereafter 
A = 0. Thus, in k-space the velocity u is given by 



tti(k) = (fc/fc_L) 2 [e m P im (k)w z (k) - i(exk)iw(k)/fc 2 



(Bl) 



Multiplying Eq. (Bl) for tii(ki) by Uj (k.2) and averaging 
over the turbulent velocity field we obtain 

4 0) (k) = (fc/fc ± ) 4 [/»/(°)(k)e mnJ R i „ i (k)P,„(k) 

+ (exk) J (exk) J G (0) (k)/fc 4 ] , (B2) 

where we assumed that the turbulent velocity field in the 
background turbulent convection is non-helical. Now we 
use an identity 



e mn Pi m (k)_P,- n (k) — eij + kij 



Py(k)-J#(kj 



(B3) 



which can be derived from 



k z (k z eij + &iki + ejk i ) 



kij k 



k^k 2 



f^(k) = uw(k)p ij (\ s .)/8wk 2 



(B5) 



In the isotropic two-dimensional turbulent flow 
G(°)(k)/fc 2 > /»/ (0) (k) and the correlation func- 
tion is given by 



^(k)=G(°)(k)P i jr(kx)/fci 



(B6) 



A simplest generalization of these correlation functions 
is an assumption that G(°)(k)/[fc 2 .f*/( )(k)] - 1 = e = 

const, and thus the correlation function /^(k) is given 



by Eq. ([121). This correlation function can be considered 
as a combination of Eqs. ( p35| ) and (B6) for the three- 
dimensional and two-dimensional turbulence. When e 
depends on the wave vector k, the correlation function 
fj^(k) is determined by two spectral functions. 

Now we derive Eq. (O) for the turbulent flux of en- 
tropy. Multiplying Eq. (|B1|) written for mQs.2) by s(ki) 
and averaging over turbulent velocity field we obtain 
Eq. (H). Multiplying Eq. (||) by i(ki_xe) 4 we get 



F(°)(k) =i(k_LXe)-*i 0) (k) 



(B7) 



Now we assume that < fr^ (k) cx <J?^/(°)(k). The integra- 



tion in k-space in Eq. flB7|) yields the numerical factor in 
Eq. (15). Note that for simplicity we assumed that the 
correlation functions F^°)(k) and /^(k) have the same 
spectrum. If these functions have different spectra, it re- 
sults only in a different value of a numerical coefficient 
in Eq. fljj). 

Now let us discuss the physical meaning of the pa- 
rameter a. To this end we derived the equation for the 
two-point correlation function $i°^(r) = (s(x) u(x + r)) 
of the turbulent flux of entropy for the background tur- 
bulent convection (which corresponds to Eq. ( |l4| ) written 
in k-space). Let us we rewrite Eq. (|lj) in the following 
form: 

$i°)(k) = (<r-e)[fc 2 +r(e.k) 2 ]<Mfc) , (B8) 
$ w (k) = ~(3-a)W(k)/8nk 4 , (B9) 

wher e T = 3 (a — l)/(3 — a). The Fourier transform of 
Eq. |H) reads 

*<°>(r) = (**-e)[A + r(e- V) 2 ]$ lu (r) , (BIO) 
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where 3> w (r) is the Fourier transform of the function 
<t w (k). Now we use the identity 



(Bll) 



where ip(r) = r 1 $J 1J (r) and ip'(r) — dip /dr. Equations 
( B10| ) and ( |B11 ) yield the two-point correlation function 



$i 0) (r) 



$(o)( r ) = ($*. e )U(r) 

,,. , 1 +rcos 2 6i 
+rip (r)- 



3 + r 



(B12) 



where 6 is the angle between e and r. The function 
ip(r) has the following properties: ip(r = 0) = 1 and 
{rip')r=o = 0, e.g., the function ip(r) = 1 — (r/Z ) 9_1 sat- 
isfies the above properties, where 1 < q < 3. Thus, the 
two-point correlation function (r) of the flux of en- 
tropy for the background turbulent convection is given 
by ' 



*(°>(i 



(g-i)(i + rcos 2 
3 + r 



9-1 



(B13) 



Simple analysis shows that — B/(q — 1) < a < 3, where 
we took into account that d&i°^ (r) /dr < for all angles 
9. The parameter a can be presented as a = [1 + £(g + 
iVCe-WCL+C/ 3 ) andC = (k/Q 9_1 -l, where l± and 
Z z are the horizontal (9 = n/2) and the vertical (9 = 0) 
scales in which the correlation function <&i°^ (r) tends to 
zero. The parameter £ describes the degree of thermal 
anisotropy. In particular, when l± = l z the parameter 
£ = and a = 1. For ij_ -C Z z the parameter £ = — 1 and 
a = — 3/(g — 1). The maximum value £ max of the param- 
eter £ is given by £ max = <Z — 1 for a = 3. Thus, for a < 1 
the thermal structures have the form of column or ther- 
mal jets (l± < lz), and a > 1 there exist the "pancake" 
thermal structures (l± > l z ) in the background turbulent 
convection. 



kij sin 2 9dCl = (87r/15)(25 ij - e y ) , 

, sin 2 6»cK7 = (87r/105)(3A ijm „ - T ijmn ) , 
/" fc^rfO = 27rPy(e) , 

kijk mn dCl — (7r/3) [-Py (e) (5 Ifm + e rml ) 
+P m (e)P mj (e) + P im (e)P nj (e)] , 
fcj kjk mn kj_ dCl = (27r/3)[Pi ra (e)e m + Pj m (e)e n ] , 

klkjk mn k~ 2 dCl = (47r/3) [P in (e)e m + P im (e)e„] , 



where Pjj(e) — (5^- e.(j, kij nin — kikjk nl k n /k^, dCl — 
sin 9 dd dip, and 



A 



ijmn 



- tjmn 



jijm 
€"nX*ijmn 



Py(k) + ei^(k ± ) = (1 + e)^ - ee, 



The above identities allow us to calculate the following 
integrals (which were used for the derivation of equations 
for the Reynolds stresses and turbulent flux of entropy): 



J r(fc)%$^(k)dk = (** • e)(r /30)[15%e* 



-10a e« 



(2a + 3)j ijm + 6bij m ] , 



r(fc)fc mn / i y(k) dk = {f*T /6)[(e/4)(T ijmn - e ijmn ) 
-(1/5 + e/4)Ai im „ + (1 + e)5ij5 

run 



APPENDIX C: THE IDENTITIES USED FOR 
THE INTEGRATION IN k-SPACE 



where 



To inte grate over the angles in k-space of Eqs. ( A3C ) 
and (A31) we used the following identities: 



kijmn dQ. = (47r/15)Aj 



<%($I)m + id-mi (**xe)j +s jm i (**xe)i]ei 



and we used an identity e q e pqn (^^ xe) m A 



ijmn — ^ijp- 
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